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A model-independent search for direct CP violation in the Cabibbo suppressed decay D + — > 
K~ K + n + in a sample of approximately 370,000 decays is carried out. The data were collected by the 
LHCb experiment in 2010 and correspond to an integrated luminosity of 35 pb _1 . The normalized 
Dalitz plot distributions for D + and D~ are compared using four different binning schemes that are 
sensitive to different manifestations of CP violation. No evidence for CP asymmetry is found. 
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I. INTRODUCTION 

To date CP violation (CPV) has been observed only 
in decays of neutral K and B mesons. All observations 
are consistent with CPV being generated by the phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [HE] of 
the Standard Model (SM). In the charm sector, CKM dy- 
namics can produce direct CP asymmetries in Cabibbo 



suppressed D ± decays of the order of 10 -3 or less [3]. 
Asymmetries of up to around 1% can be generated by 
new physics (NP) [U [5]. In most extensions of the SM, 
asymmetries arise in processes with loop diagrams. How- 
ever, in some cases CPV could occur even at tree level, 
for example in models with charged Higgs exchange. 

In decays of hadrons, CPV can be observed when 
two different amplitudes with non-zero relative weak 
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and strong phases contribute coherently to a final state. 
Three-body decays are dominated by intermediate reso- 
nant states, and the requirement of a non-zero relative 
strong phase is fulfilled by the phases of the resonances. 
In two-body decays, CPV leads to an asymmetry in the 
partial widths. In three-body decays, the interference 
between resonances in the two-dimensional phase space 
can lead to observable asymmetries which vary across the 
Dalitz plot. 

CP- violating phase differences of 10° or less do not, 
in general, lead to large asymmetries in integrated decay 
rates, but they could have clear signatures in the Dalitz 
plot, as we will show in Sect. |III| This means that a two- 
dimensional search should have higher sensitivity than an 
integrated measurement. In addition, the distribution 
of an asymmetry across phase space could hint at the 
underlying dynamics. 

At present, no theoretical tools for computing decay 
fractions and relative phases of resonant modes in D de- 
cays have been applied to multibody D + decay modes, 
and no predictions have been made for how asymmetries 
might vary across their Dalitz plots. A full Dalitz plot 
analysis of large data samples could, in principle, mea- 
sure small phase differences. However, rigorous control 
of the much larger strong phases would be required. For 
this to be achieved, better understanding of the ampli- 
tudes, especially in the scalar sector, will be needed, and 
effects like three-body final state interactions should be 
taken into account. 

This paper describes a model-independent search for 
direct CPV in the Cabibbo suppressed decay D + — > 
K~ K + tt + in a binned Dalitz plot[J A direct comparison 
between the D + and the D~ Dalitz plots is made on a 
bin-by-bin basis. The data sample used is approximately 
35 pb _1 collected in 2010 by the LHCb experiment at a 
centre of mass energy of y/s = 7 TeV. This data set corre- 
sponds to nearly 10 and 20 times more signal events than 
used in previous studies of this channel performed by the 
BABAR [H] and CLEO-c [7] collaborations, respectively. 
It is comparable to the dataset used in a more recent 
search for CPV in D+ -> <f>ir+ decays at BELLE 0. 

The strategy is as follows. For each bin in the Dalitz 
plot, a local CP asymmetry variable is defined [51 IIP). 

s% = N l {D+)-aN\D-) = N tot (D+) 

CP y/N*(D+) + a 2 N*{D-) ' a N tot (D-)' 

(!) 

where N l (D + ) and N l (D ) are the numbers of D ± can- 
didates in the i th bin and a is the ratio between the 
total D + and D~ yields. The parameter a accounts for 
global asymmetries, i.e. those that are constant across 
the Dalitz plot. 



1 Throughout this paper charge conjugation is implied, unless oth- 
erwise stated. 



In the absence of Dalitz plot dependent asymmetries, 
the S l CP values are distributed according to a Gaussian 
distribution with zero mean and unit width. CPV signals 
are, therefore, deviations from this behaviour. The nu- 
merical comparison between the D + and the D~ Dalitz 
plots is made with a x 2 test (x 2 = XX^cp) 2 )- The num- 
ber of degrees of freedom is the number of bins minus one 
(due to the constraint of the overall D + /D~ normaliza- 
tion). The p- value that results from this test is defined 
as the probability of obtaining, for a given number of de- 
grees of freedom and under the assumption of no CPV, 
a x 2 that is at least as high as the value observed [llj . 
It measures the degree to which we are confident that 
the differences between the D + and D~ Dalitz plots are 
driven only by statistical fluctuations. 

If CPV is observed, the p- value from this test could be 
converted into a significance for a signal using Gaussian 
statistics. However, in the event that no CPV is found, 
there is no model-independent mechanism for setting lim- 
its on CPV within this procedure. In this case, the results 
can be compared to simulation studies in which an arti- 
ficial CP asymmetry is introduced into an assumed am- 
plitude model for the decay. Since such simulations are 
clearly model-dependent, they are only used as a guide 
to the sensitivity of the method, and not in the determi- 
nation of the p-values that constitute the results of the 
analysis. 

The technique relies on careful accounting for local 
asymmetries that could be induced by sources such as, 
the difference in the if-nucleon inelastic cross-section, 
differences in the reconstruction or trigger efficiencies, 
left-right detector asymmetries, etc. These effects arc in- 
vestigated in the two control channels D + — > K~tt + tt + 
and D+ -> K-K+TT+. 

The optimum sensitivity is obtained with bins of nearly 
the same size as the area over which the asymmetry ex- 
tends in the Dalitz plot. Since this is a search for new 
and therefore unknown phenomena, it is necessary to be 
sensitive to effects restricted to small areas as well as 
those that extend over a large region of the Dalitz plot. 
Therefore two types of binning scheme are employed. 
The first type is simply a uniform grid of equally sized 
bins. The second type takes into account the fact that the 
D + — > K~K + ir + Dalitz plot is dominated by the <f>ir + 
and K*(892)°K + resonances, so the event distribution 
is highly non-uniform. This "adaptive binning" scheme 
uses smaller bins where the density of events is high, aim- 
ing for a uniform bin population. In each scheme, differ- 
ent numbers of bins are used in our search for localized 
asymmetries. 

The paper is organized as follows. In Sect.|Hj a descrip- 
tion of the LHCb experiment and of the data selection 
is presented. In Sect. |III| the methods and the binnings 
are discussed in detail. The study of the control channels 
and of possible asymmetries generated by detector effects 
or backgrounds is presented in Sect. |IV[ The results of 
our search are given in Sect. [Vj and the conclusions in 
Sect. EI 
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II. DETECTOR, DATASET AND SELECTION 

The LHCb detector [12] is a single-arm forward spec- 
trometer with the main purpose of measuring CPV and 
rare decays of hadrons containing b and c quarks. A ver- 
tex locator (VELO) determines with high precision the 
positions of the vertices of primary pp collisions (PVs) 
and the decay vertices of long-lived particles. The track- 
ing system also includes a large area silicon strip detec- 
tor located in front of a dipole magnet with an inte- 
grated field of around 4 Tm, and a combination of silicon 
strip detectors and straw drift chambers placed behind 
the magnet. Charged hadron identification is achieved 
with two ring-imaging Cherenkov (RICH) detectors. The 
calorimeter system consists of a preshower, a scintilla- 
tor pad detector, an electromagnetic calorimeter and a 
hadronic calorimeter. It identifies high transverse en- 
ergy (Et) hadron, electron and photon candidates and 
provides information for the trigger. Five muon sta- 
tions composed of multi-wire proportional chambers and 
triple-GEMs (gas electron multipliers) provide fast infor- 
mation for the trigger and muon identification capability. 

The LHCb trigger consists of two levels. The first, 
hardware-based level selects leptonic and hadronic final 
states with high tranverse momentum, using the subset 
of the detectors that are able to reduce the rate at which 
the whole detector is read out to a maximum of 1 MHz. 
The second level, the High Level Trigger (HLT), is sub- 
divided into two software stages that can use the infor- 
mation from all parts of the detector. The first stage, 
HLT1, performs a partial reconstruction of the event, 
reducing the rate further and allowing the next stage, 
HLT2, to fully reconstruct the individual channels. At 
each stage, several selections designed for specific types 
of decay exist. As luminosity increased throughout 2010 
several changes in the trigger were required. To match 
these, the datasets for signal and control modes are di- 
vided into three parts according to the trigger, samples 
1, 2 and 3, which correspond to integrated luminosities 
of approximately 3, 5 and 28 pb _1 respectively. The 
magnet polarity was changed several times during data 
taking. 

The majority of the signal decays come via the 
hadronic hardware trigger, which has an Et threshold 
that varied between 2.6 and 3.6 GeV in 2010. In the 
HLT1, most candidates also come from the hadronic se- 
lections which retain events with at least one high trans- 
verse momentum (pt) track that is displaced from the 
PV. In the HLT2, dedicated charm triggers select most 
of the signal. However, the signal yield for these channels 
can be increased by using other trigger selections, such as 
those for decays of the form B — > DX. To maintain the 
necessary control of Dalitz plot-dependent asymmetries, 
only events from selections which have been measured 
not to introduce charge asymmetries into the Dalitz plot 
of the D + — » K~ir + Tr + control mode are accepted. 

The signal ( D+ -> R-K+TT+) and control (D+ -> 
K~7r + ir + and —> K~ K + n + ) mode candidates are 



selected using the same criteria, which are chosen to max- 
imize the statistical significance of the signal. Moreover, 
care is taken to use selection cuts that do not have a low 
efficiency in any part of the Dalitz plot, as this would re- 
duce the sensitivity in these areas. The selection criteria 
are the same regardless of the trigger conditions. 

The event selection starts by requiring at least one PV 
with a minimum of five charged tracks to exist. To con- 
trol CPU consumption each event must also have fewer 
than 350 reconstructed tracks. The particle identification 
system constructs a relative log-likelihood for pion and 
kaon hypotheses, DLLr-^, and we require DLLa'tt > 7 
for kaons and < 3 for pions. Three particles with appro- 
priate charges are combined to form the candidates. 
The corresponding tracks are required to have a good 
fit quality (x 2 /ndf < 5), pr > 250 MeV/c, momentum 
p > 2000 MeV/c and the scalar sum of their pr above 
2800 MeV/c. Because a typical D + travels for around 
8 mm before decaying, the final state tracks should not 
point to the PV. The smallest displacement from each 
track to the PV is computed, and a x 2 (Xip)j formed by 
using the hypothesis that this distance is equal to zero, 
is required to be greater than 4 for each track. The three 
daughters should be produced at a common origin, the 
charm decay vertex, with vertex fit x 2 / n df < 10. 

This 'secondary' vertex must be well separated from 
any PV, thus a flight distance variable (xfd) ^ s con " 
structed. The secondary vertex is required to have 
x| D > 100, and to be downstream of the PV. The pt 
of the Df, candidate must be greater than 1000 MeV/c, 
and its reconstructed trajectory is required to originate 
from the PV (x 2 P < 12). 

In order to quantify the signal yields (S), a simultane- 
ous fit to the invariant mass distribution of the D + and 
D~ samples is performed. A double Gaussian is used for 
the A^~A' + 7r + signal, whilst the background (B) is de- 
scribed by a quadratic component and a single Gaussian 
for the small contamination from D* + — > D°(K~ K + )ir + 
above the D+ peak. The fitted mass spectrum for sam- 
ples 1 and 3 combined is shown in Fig. [T] giving the 
yields shown in Table [I] A weighted mean of the widths 
of the two Gaussian contributions to the mass peaks is 
used to determine the overall widths, a, as 6.35 MeV/c 2 
for D+ -> K-K+1T+, 7.05 MeV/c 2 for D+ -> K-R+TT+, 
and 8.0 MeV/c 2 for D + — > K~tt + tt + . These values are 
used to define signal mass windows of approximately 2a 
in which the Dalitz plots are constructed. The purities, 
defined as S/ (B + S) within these mass regions, are also 
shown in Table [T] for samples 1 and 3 in the different 
decay modes. 

For sample 2, the yield cannot be taken directly from 
the fit, because there is a mass cut in the HLT2 line that 
accepts the majority of the signal, selecting events in a 
±25 MeV/c 2 window around the nominal value. How- 
ever, another HLT2 line with a looser mass cut that 
is otherwise identical to the main HLT2 line exists, al- 
though only one event in 100 is retained. In this line the 
purity is found to be the same in sample 2 as in sample 
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FIG. 1. Fitted mass spectra of (a) K tt + tt + and (b) K K + ir + candidates from samples 1 and 3, D + and D combined. The 
signal mass windows and sidebands defined in the text are labelled. 
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FIG. 2. Dalitz plot of the D + -> K~K + tv + decay for se- 
lected candidates in the signal window. The vertical K* (892)° 
and horizontal (^(1020) contributions are clearly visible in the 
data. 



3. The yield in sample 2 is then inferred as the total 
(S + B) in all allowed triggers in the mass window times 
the purity in sample 3. Thus the overall yield of signal 
D + — > K~ K + tt + candidates in the three samples within 
the mass window is approximately 370,000. The total 
number of candidates (S + B) in each decay mode used 
in the analysis are given in Table [TTJ The Dalitz plot of 
data in the D + window is shown in Fig. [2j 

Within the 2a D + — > K~K + tt + mass window, about 
8.6% of events are background. Apart from random 
three-body track combinations, charm backgrounds and 
two-body resonances plus one track are expected. Charm 
reflections appear when a particle is wrongly identified in 
a true charm three-body decay and/or a track in a four- 
body charm decay is lost. The main three-body reflec- 
tion in the K~ K + n + spectrum is the Cabibbo- favoured 



TABLE I. Yield (5) and purity for samples I and 3 after 
the final selection. The purity is estimated in the 2a mass 
window. 





Decay 


Yield 


Purity 






Sample 1+3 


Sample 1 Sample 3 


D+ 


-> K-K+TT+ 


(3.284 ± 0.006) x 1 b 


88% 92% 


Dt 


-> K~K + iT + 


(4.615 ±0.012) x 10 5 


89% 92% 


D+ 




(3.3777 ± 0.0037) x 10 6 


98% 98% 



TABLE II. Number of candidates (S + B) in the signal win- 
dows shown in Fig. [I] after the final selection, for use in the 
subsequent analysis. 

sample 1 sample 2 sample 3 Total 



D+- 
D+ 



K~K + TT + 84,667 65,781 253,446 403,894 
K-K+-K+ 126,206 91,664 346,068 563,938 
• R--K+-K+ 858,356 687,197 2,294,315 3,839,868 



D + — > K~tt + it + , where the incorrect assignment of the 
kaon mass to the pion leads to a distribution that par- 
tially overlaps with the — > K~K + ir + signal region, 
but not with D + — > K~ K + n + . The four body, Cabibbo- 
favoured mode D° —> K~ir + ir~ir + where a ir + is lost 
and the n~ is misidcntificd as a K~ will appear broadly 
distributed in K~K + n + mass, but its resonances could 
create structures in the Dalitz plot. Similarly, X*(892)° 
and <f) resonances from the PV misreconstructed with a 
random track forming a three-body vertex will also ap- 
pear. 
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TABLE III. The CLEO-c amplitude model "B" [7j used in the simulation studies. The uncertainties are statistical, experimental 
systematic and model systematic respectively. 



Resonance 



Amplitude 



Relative phase 



A"* (892)° 

F5(1430)° 

/<800) 

K* 2 (U30)° 

0(1020) 

a (1450)° 

0(1680) 



Fit fraction 

r+0. 4+0.1 



1 (fixed) 

+ n 1Q+0. 10+0.42 
4.00 =t U.li_ 01 _ Q 39 

2.30±0.13±g-?i±S-g 
7.6±0.8t^ 

1.166 ±0.015™t°o:°o 2 o 5 9 
1.50±0.10± 09+a92 



1.86 ±0.20 



0.06-0.33 
+0.02+0.62 
-0.08-0.77 



(fixed) 
70 ±6tjt*» 
-87±6l 2 3 li 5 
171 ± 

-163 ±3t\tl i 
116 ± 2±1±^ 



25.7 iO.5" 1 

18.8 ± 1.2±° : f±| : | 

7.o±o.8±8:g±?;s 

l.7±0.4±^± ; 2 7 
27.8 ±0.412™ 
4.6±0.6±°-§jl™ 

51 +0 1 T +0.01+0.3T 
U.Ol =t= U-ll_o.Q4_o.i5 



III. METHODS AND BINNINGS 

Monte Carlo pseudo-experiments are used to verify 
that we can detect CPV with the strategy outlined in 
Sect, [ijwithout producing fake signals, and to devise and 
test suitable binning schemes for the Dalitz plot. They 
are also used to quantify our sensitivity to possible man- 
ifestations of CPV, where we define the sensitivity to a 
given level of CPV as the probability of observing it with 
3(7 significance. 

For the D + — > K~K + ir + Dalitz plot model, the result 
of the CLEO-c analysis (fit B) [7] is used. The ampli- 
tudes and phases of the resonances used in this model 
are reproduced in Table |LTl} For simplicity, only resonant 
modes with fit fractions greater than 2% are included 
in the pseudo-experiments. The fit fraction for a reso- 
nance is defined as the integral of its squared amplitude 
over the Dalitz plot divided by the integral of the square 
of the overall complex amplitude. An efficiency func- 
tion is determined from a two-dimensional second order 
polynomial fit to the Dalitz plot distribution of triggered 
events that survive the selection cuts in the GEANT- 
based |T3] LHCb Monte Carlo simulation for nonresonant 
D + — > K~K + tt + . A simple model for the background 
is inferred from the Dalitz plots of the sidebands of the 
D + — > K~ K + ir + signal. It is composed of random com- 
binations of K~ , K + , and n + tracks, <j) resonances with 
7T + tracks, and K* (892)° resonances with K + tracks. 
The CLEO-c Dalitz plot analysis has large uncertainties, 
as do the background and efficiency simulations (due to 
limited numbers of MC events) , so the method is tested 
on a range of different Dalitz plot models. 

Pseudo-experiments with large numbers of events are 
used to investigate how CPV would be observed in the 
Dalitz plot. These experiments are simple "toy" simu- 
lations that produce points in the Dalitz plot according 
to the probability density function determined from the 
CLEO-c amplitude model with no representation of the 
proton-proton collision, detector, or trigger. Figure [3^a) 
illustrates the values of S l CP observed with 8 x 10 7 events 
and no CPV. This dataset is approximately 50 times 
larger than the data sample under study. The result- 
ing x 2 /ndf is 253.4/218, giving a p-value for consistency 
with no CPV of 5.0%. This test shows that the method 



by itself is very unlikely to yield false positive results. 
Figure J3^b) shows an example test using 5 x 10 7 events 
with a CP violating phase difference of 4° between the 
amplitudes for the 0(lO2O)7r + component in D + and D~ 
decays. The p-valuc in this case is less than 10~ lo °. The 
CPV effect is clearly visible, and is spread over a broad 
area of the plot, changing sign from left to right. This 
sign change means the CPV causes only a 0.1% differ- 
ence in the total decay rate between D + and D~ . This 
illustrates the strength of our method, as the asymme- 
try would be much more difficult to detect in a measure- 
ment that was integrated over the Dalitz plot. Even with 
no systematic uncertainties, to see a 0.1% asymmetry at 
the 3a level would require 2.25 x 10 6 events. With the 
method and much smaller dataset used here we would 
observe this signal at the 3cr level with 76% probability, 
as shown in Table ITVl below. 

The sensitivity to a particular manifestation of CPV 
depends on the choice of binning. The fact that the CP- 
violating region in most of the pseudo-experiments covers 
a broad area of the Dalitz plot suggests that the optimal 
number of bins for this type of asymmetry is low. Each 
bin adds a degree of freedom without changing the x 2 
value for consistency with no CPV. However, if CP asym- 
metries change sign within a bin, they will not be seen. 
Similarly, the sensitivity is reduced if only a small part 
of a large bin has any CPV in it. To avoid effects due 
to excessive fluctuations, bins that contain fewer than 50 
candidates are not used anywhere in the analysis. Such 
bins are very rare. 

The binnings are chosen to reflect the highly non- 
uniform structure of the Dalitz plot. A simple adap- 
tive binning algorithm was devised to define binnings of 
approximately equal population without separating D + 
and D~ . Two binnings that are found to have good 
sensitivity to the simulated asymmetries contain 25 bins 
( "Adaptive I" ) arranged as shown in Fig. Qa) , and 106 
bins ("Adaptive II") arranged as shown in Fig.|4^b). For 
Adaptive I, a simulation of the relative value of the strong 
phase across the Dalitz plot in the CLEO-c amplitude 
model is used to refine the results of the algorithm: if 
the strong phase varies significantly across a bin, CP 
asymmetries are more likely to change sign. Therefore 
the bin boundaries are adjusted to minimise changes in 
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FIG. 3. Sep across the Dalitz plot in a Monte Carlo pseudo-experiment with a large number of events with (a) no CPV and 
(b) a 4° CPV in the 4>tt phase. Note the difference in colour scale between (a) and (b). 
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FIG. 4. Layout of the (a) "Adaptive I" and (b) "Adaptive II" binnings on the Dalitz plot of data. 



the strong phase within bins. The model-dependence of 
this simulation could, in principle, influence the binning 
and therefore the sensitivity to CPV, but it cannot intro- 
duce model-dependence into the final results as no artifi- 
cial signal could result purely from the choice of binning. 
Two further binning schemes, "Uniform I" and "Uniform 
IF , are defined. These use regular arrays of rectangular 
bins of equal size. 

The adaptive binnings are used to determine the sen- 
sitivity to several manifestations of CPV. With 200 test 
experiments of approximately the same size as the sig- 
nal sample in data, including no asymmetries, no CP- 
violating signals are observed at the 3<r level with Adap- 
tive I or Adaptive II. The expectation is 0.3. 

With the chosen binnings, a number of sets of 100 
pseudo-experiments with different CP- violating asymme- 
tries are produced. The probability of observing a given 
signal in either the 0(1020) or k(800) resonances with 
3cr significance is calculated in samples of the same size 
as the dataset. The results are given in Table |IV| The 
CPV shows up both in the x 2 / ncu? and in the width of 



the fitted Sep distribution. 

For comparison, the asymmetries in the <f> phase and 
k magnitude measured by the CLEO collaboration us- 
ing the same amplitude model were (6 ± 6^2-2)° anc ^ 
( — 12 ± 12j^! 10 )% 2 where the uncertainties are statisti- 
cal, systematic ana model-dependent, respectively. Ta- 
ble [TV] suggests that, assuming their model, we would be 
at least 95% confident of detecting the central values of 
these asymmetries. 

The sensitivity of the results to variations in the Dalitz 
plot model and the background is investigated, and ex- 
ample results for the CP asymmetry in the 0(1020) phase 
are shown in Table [V] In this table, models A and B 
are taken from the CLEO paper, model B2 includes an 
/o(980) contribution that accounts for approximately 8% 



2 The conventions used in the CLEO paper to define asymmetry 
are different, so the asymmetries in Table II of [7] have been 
multiplied by two in order to be comparable with those given 
above. 
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TABLE IV. Results from sets of 100 pseudo-experiments with different CP asymmetries and Adaptive I and II binnings. p(3a) 
is the probability of a 3a observation of CPV. (S) is the mean significance with which CPV is observed. 



CPV Adaptive I Adaptive II 





p(3a) 


(S) 


P(3a) 


(S) 


no CPV 





0.84a 


1% 


0.84cr 


6° in 0(1020) phase 


99% 


7.0a 


98% 


5.2a 


5° in 0(1020) phase 


97% 


5.5a 


79% 


3.8a 


4° in 0(1020) phase 


76% 


3.8a 


41% 


2.7a 


3° in 0(1020) phase 


38% 


2.8a 


12% 


1.9a 


2° in 0(1020) phase 


5% 


1.6a 


2% 


1.2a 


6.3% in k(800) magnitude 


16% 


1.9a 


24% 


2.2a 


11% in k(800) magnitude 


83% 


4.2a 


95% 


5.6a 



TABLE V. Results from sets of 100 pseudo-experiments with 4° CPV in the 0(1020) phase and different Dalitz plot models. 
p(3a) is the probability of a 3a observation of CPV. (S) is the mean significance with which CPV is observed. The sample size 
is comparable to that seen in data. 



Model Adaptive I Adaptive II 





p(3a) 


(S) 


p(3a) 


(S) 


B (baseline) 


76% 


3.8a 


41% 


2.7a 


A 


84% 


4.3a 


47% 


2.9a 


B2 (add /o(980)) 


53% 


3.2a 


24% 


2.2a 


B3 (vary j?S(1430) magn.) 


82% 


4.0a 


41% 


2.8a 


B4 (vary i?S(1430) phase) 


73% 


3.7a 


38% 


2.7a 



of events, and models B3 and B4 are variations of the 
Kq(1A30)° amplitude and phase within their uncertain- 
ties. As expected, the sensitivity to CPV in the reso- 
nances of an amplitude model depends quite strongly on 
the details of the model. This provides further justifi- 
cation for our model-independent approach. However, a 
reasonable level of sensitivity is retained in all the cases 
we tested. Thus, when taken together, the studies show 
that the method works well. It does not yield fake signals, 
and should be sensitive to any large CPV that varies sig- 
nificantly across the Dalitz plot even if it does not occur 
precisely in the way investigated here. 



IV. CONTROL MODES 

It is possible that asymmetries exist in the data that 
do not result from CPV, for example due to production, 
backgrounds, instrumental effects such as left-right dif- 
ferences in detection efficiency, or momentum-dependent 
differences in the interaction cross-sections of the daugh- 
ter particles with detector material. Our sensitivity to 
such asymmetries is investigated in the two Cabibbo 
favoured control channels, where there is no large CPV 
predicted. The D + — » K~ii + ii + control mode has an 
order of magnitude more candidates than the Cabibbo- 
suppressed signal mode, and is more sensitive to detec- 
tor effects since there is no cancellation between K + and 
K~ reconstruction efficiencies. Conversely, the Df —> 



K~ K +r K + control mode is very similar to our signal mode 
in terms of resonant structure, number of candidates, 
kinematics, detector effects, and backgrounds. 

The control modes and their mass sidebands defined 
in Fig. [I] are tested for asymmetries using the method 
described in the previous section. Adaptive and uniform 
binning schemes are denned for D + — > K~-k + ii + and 
Df — > K~K + tt + . They are applied to samples 1-3 and 
each magnet polarity separately. In the final results, the 
asymmetries measured in data taken with positive and 
negative magnet polarity are combined in order to cancel 
left-right detector asymmetries. The precise number of 
bins chosen is arbitrary, but care is taken to use a wide 
range of tests with binnings that reflect the size of the 
dataset for the decay mode under study. 

For D + — > K~tt + it + , five different sets of bins in each 
scheme are used. A very low p-value would indicate a 
local asymmetry. One test with 25 adaptive bins in one 
of the subsamples (with negative magnet polarity) has 
a p-value of 0.1%, but when combined with the positive 
polarity sample the p-value increases to 1.7%. All other 
tests yield p- values ranging from 1-98%. Some example 
results are given in Table |VI| A typical distribution of 
the Sqp values with a Gaussian fit is shown in Fig. [5^ a) 
for a test with 900 uniform bins. The fitted values of 
the mean and width are consistent with one and zero 
respectively, suggesting that the differences between the 
D + and the D~ Dalitz plots are driven only by statistical 
fluctuations. 
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FIG. 5. (a) Distribution of Sep values from D + — > K tt + -k + from a test with 900 uniform bins. The mean of the fitted 
Gaussian distribution is 0.015 ± 0.034 and the width is 0.996 ± 0.023. (b) Distribution of Sep values from D+ -> K~K + tt + 
with 129 bins. The fitted mean is -0.011 ± 0.084 and the width is 0.958 ± 0.060. 
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FIG. 6. Dalitz plots of (a) D + — > K~tx + -k + , showing the 25-bin adaptive scheme with the Sep values, and (b) Df — > K~ K + ix + , 
showing the three regions referred to in the text. The higher and lower K~n + invariant mass combinations are plotted in (a) 
as there are identical pions in the final state. 



TABLE VI. Results (p-values, in %) from tests with the 
D + —¥ K~~-k + tt + control channel using the uniform and adap- 
tive binning schemes. The values correspond to tests per- 
formed on the whole dataset in the mass windows defined in 
Sect. HIl 



1300 bins 900 bins 400 bins 100 bins 25 bins 
Uniform 7X8 1777 TU> 54l3 1.7 
Adaptive 81.7 57.4 65.8 30.0 11.8 



TABLE VII. Results (p-values, in %) from tests with the 
Df — > K~ K + tt + control channel using the uniform binning 
scheme. The values correspond to tests performed separately 
on Zones A-C, with samples 1-3 and both magnet polarities 
combined. 



bins Zone A Zone B Zone C 



300 


20.1 


25.3 


14.5 


100 


41.7 


84.6 


89.5 


30 


66.0 


62.5 


24.6 



For the Df K~ K + ir + mode a different procedure 
is followed due to the smaller sample size and to the high 
density of events along the </> and the K*(892)° bands. 
The Dalitz plot is divided into three zones, as shown in 
Fig. [6] Each zone is further divided into 300, 100 and 30 
bins of same size. The results are given in Table |VII| In 
addition, a test is performed on the whole Dalitz plot us- 



ing 129 bins chosen by the adaptive algorithm, and a ver- 
sion of the 25-bin scheme outlined in Sect. |III| scaled by 
the ratio of the available phase space in the two modes. 
These tests yield p-values of 71.5% and 34.3% respec- 
tively. 

Other possible sources of local charge asymmetry in the 
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TABLE VIII. Results (p-values, in %) from tests with the 
D + — > K~ir + ir + and Df — > K~K + -k + samples divided ac- 
cording to the impact parameter with respect to the primary 
vertex. The tests are performed using the adaptive binning 
scheme with 25 bins. 



CPV is investigated with uniform binnings, using 200 
and 530 bins of equal size. For each of these binning 
choices, the significance S l CP of the difference in D + and 







xip < 6 


xip >6 


D+ 




8.5 


88.9 


Df- 


-> K-K+TY+ 


52.0 


30.6 



signal region are the charm contamination of the back- 
ground, and asymmetries from CPV in misreconstructed 
B decays. In order to investigate the first possibility, 
similar tests are carried out in the mass sidebands of the 
— > K~K + n + signal (illustrated in Fig.jlJ). There is 
no evidence for asymmetries in the background. 

From a simulation of the decay D + — s- K~tt + it + the 
level of secondary charm (B — > DX) in our selected sam- 
ple is found to be 4.5%. The main discriminating vari- 
able to distinguish between prompt and secondary charm 
is the impact parameter (IP) of the D with respect to the 
primary vertex. Given the long B lifetime, the IP distri- 
bution of secondary charm candidates is shifted towards 
larger values compared to that of prompt D + mesons. 

The effect of secondary charm is investigated by divid- 
ing the data set according to the value of the candidate 
IP significance (xjp)- The subsample with events having 
larger xjp are likely to be richer in secondary charm. The 
results are shown in Table IVlTll No anomalous effects are 
seen in the high xjp sample, so contamination from sec- 
ondary charm with CPV does not affect our results for 
studies with our current level of sensitivity. 

The analysis on the two control modes and on the 
sidebands in the final states K~K + ir + and K~ir + Tr + 
gives results from all tests that are fully consistent with 
no asymmetry. Therefore, any asymmetry observed in 
D + — > K~ K + ir + is likely to be a real physics effect. 



V. RESULTS 

The signal sample with which we search for CP vio- 
lation consists of 403,894 candidates selected within the 
K-K+-K+ mass window from 1856.7 to 1882.1 MeV/c 2 , 
as described in Sect. [Il] There are 200,336 and 203,558 
D + and D~ candidates respectively. This implies a nor- 
malization factor a = N tot (D + )/N tot (D^) = 0.984 ± 
0.003, to be used in Eq. [T] 

The strategy for looking for signs of localized CPV is 
discussed in the previous sections. In the absence of local 
asymmetries in the control channels D + K~tt + tt + and 
£)+ — > K~K + ir + and in the sidebands of the K~ K + iv + 
mass spectrum, we investigate the signal sample under 
different binning choices. 

First, the adaptive binning is used with 25 and 106 
bins in the Dalitz plot as illustrated in Fig. [4] Then 



D p< 
Eq.0 



opulation is computed for each bin i, as defined in 
The x 2 / n( ff = ^2i(Scp) 2 /ndf is calculated and 



the p-value is obtained. The distributions of S^p are 
fitted to Gaussian functions. 

The p- values are shown in Table |IX| The Dalitz plot 
distributions of S l CP are shown in Fig.[7| In Fig.[8]the dis- 
tributions of Sq P and the corresponding Gaussian fits for 
the different binnings are shown. The p-values obtained 
indicate no evidence for CPV. This is corroborated by 
the good fits of the S l CP distributions to Gaussians, with 
means and widths consistent with and 1, respectively. 

As further checks, many other binnings are tested. 
The number of bins in the adaptive and uniform bin- 
ning schemes is varied from 28 to 106 and from 21 to 
530 respectively. The samples are separated according 
to the magnet polarity and the same studies are re- 
peated. In all cases the p-values are consistent with no 
CPV, with values ranging from 4% to 99%. We conclude 
that there is no evidence for CPV in our data sample of 
D+ -> K-K + n + . 



VI. CONCLUSION 

Due to the rich structure of their Dalitz plots, three 
body charm decays are sensitive to CP violating phases 
within and beyond the Standard Model. Here, a model- 
independent search for direct CP violation is performed 
in the Cabibbo suppressed decay D + — > K~ K + n + with 
35 pb _1 of data collected by the LHCb experiment, and 
no evidence for CPV is found. Several binnings are used 
to compare normalised D + and D~ Dalitz plot distri- 
butions. This technique is validated with large num- 
bers of simulated pseudo-experiments and with Cabibbo 
favoured control channels from the data: no false positive 
signals are seen. To our knowledge this is the first time 
a search for CPV is performed using adaptive bins which 
reflect the structure of the Dalitz plot. 

Monte Carlo simulations illustrate that large localised 
asymmetries can occur without causing detectable differ- 
ences in integrated decay rates. The technique used here 
is shown to be sensitive to such asymmetries. Assuming 
the decay model, efficiency parameterisation and back- 
ground model described in Sect. Ill we would be 90% 



confident of seeing a CP violating difference of either 5° 
in the phase of the (f)ir + or 11% in the magnitude of the 
k(800)AT + with 3cr significance. Since we find no evidence 
of CPV, effects of this size are unlikely to exist. 
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TABLE IX. Fitted means and widths, x 2 /ndi and p-values for consistency with no CPV for the D + — > K K + ty + decay mode 
with four different binnings. 



Binning Fitted mean Fitted width \ /ndf p- value (%) 
Adaptive I 0.01 ± 0.23 1.13 ± 0.16 32.0/24 1277 
Adaptive II -0.024 ± 0.010 1.078 ± 0.074 123.4/105 10.6 
Uniform I -0.043 ± 0.073 0.929 ± 0.051 191.3/198 82.1 
Uniform II -0.039 ± 0.045 1.011 ± 0.034 519.5/529 60.5 
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FIG. 7. Distribution of S^ P in the Dalitz plot for (a) "Adaptive I", (b) "Adaptive II", (c) "Uniform I" and (d) "Uniform II". 
In (c) and (d) bins at the edges are not shown if the number of entries is not above a threshold of 50 (see Sect. Ill I. 
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